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Abstract 



Recent work on Lagrangian descriptors has shown that Lyapunov Exponents 
can be apphed to observed or simulated data to characterize the horizontal stirring 
and transport properties of the oceanic flow. However, a more detailed analysis 
of regional dependence and seasonal variability was still lacking. In this paper, 
we analyze the near-surface velocity field obtained from the Ocean general circula- 
tion model For the Earth Simulator (OFES) using Finite-Size Lyapunov Exponents 
(FSLE). We have characterized regional and seasonal variability. Our results show 
that horizontal stirring, as measured by FSLEs, is seasonally-varying, with maxi- 
mum values in Summer time. FSLEs also strongly vary depending on the region: 
we have first characterized the stirring properties of Northern and Southern Hemi- 
spheres, then the main oceanic basins and currents. We have finally studied the 
relation between averages of FSLE and some Eulerian descriptors such as Eddy 
Kinetic Energy (EKE) and vorticity (uj) over the different regions. 

1 Introduction 

A detailed knowledge of the transport, dispersion, stirring and mixing mechanisms of 
water masses across the global ocean is of crucial interest to fully understand, for example, 
heat and tracer budgets, or the role of oceans in climate regulation. There has been a 
recent strong activity in the study of these processes from a Lagrangian perspective. Some 
works have addressed the global variability of them using finite-time Lyapunov exponents 
(FTLEs) computed from currents derived from satellite altimetry [llll2]. These studies 
quantify stirring intensity, and identify mesoscale eddies and other Lagrangian Coherent 
Structures (LCSs). Furthermore, previous works |13] pointed out relationships between 
Lagrangian and Eulerian quantifiers of stirring/mixing activity (FTLEs and Eddy Kinetic 
Energy (EKE) or mean strain rate). 

Having in mind the implications for the distribution of biogeochemical tracers, our 
goal is to extend the previous works to provide detailed seasonal analysis and a com- 
parative study between different ocean regions and different scales: Earth's hemispheres, 
ocean basins, and boundary currents. To this end we use finite-size Lyapunov exponents 
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(FSLEs). These quantities are related to FTLEs since they also compute stretching and 
contraction time scales for transport, but they depend on exphcit spatial scales which are 
simple to specify and to interpret in oceanographic contexts flUi [TT| 120 1 HOj. In particular 
we will focus on the impact on transport of mesoscale processes, for which characteristic 
spatial scales as a function of latitude are well known. We are also interested in checking 
the existence of relationships between Lagrangian measures of horizontal stirring inten- 
sity, as given by averages of finite-size Lyapunov exponents (FSLE), and other dynamic, 
Eulerian quantities, such as EKE or vorticity. Such a functional relation does not need 
to hold in general, but may be present when there is a connection between the mecha- 
nisms giving rise to mesoscale turbulence (probably, baroclinic instability) and horizontal 
stirring. 

The paper is organized as follows. In Section [2] we describe the data and tools used 
in this study. In section |3] we first present the geographical and seasonal characterization 
of the horizontal stirring, and then we investigate the relation of FSLE with EKE and 
vorticity. Finally, in the Conclusions we present a summary and concluding remarks. 

2 Data and Methods 

Our dataset consists of an output from the Ocean general circulation model For the Earth 
Simulator (OFES) |26l [25]. This is a near-global ocean model that has been spun up 
for 50 years under climatological forcing taken from monthly mean NCEP (United States 
National Centers for Environmental Prediction) atmospheric data. After that period the 
OFES is forced by the daily mean NCEP reanalysis for 48 years from 1950 to 1998. See 
[26] for additional details on the forcing. The output of the model corresponds to daily 
data for the last 8 years. Horizontal angular resolution is the same in both the zonal, 0, 
and meridional, 6, directions, with values of = A0 = 1/10°. The output has been 
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interpolated to 54 vertical z-layers and has a temporal resolution of one day. The velocity 
fields that we have used in this work correspond to the first two years, 1990 and 1991, of 
the output. Vertical displacements are unimportant during the time scales we consider 
here so that, despite horizontal layers are not true isopycnals, most fiuid elements remain 
in their initial horizontal layer during the time of our Lagrangian computation. Thus we 
use in our analysis horizontal velocities in single horizontal layers. We refer to recent works 
[28t E] for Lyapunov analyses considering vertical displacements. Unless explicitly stated, 
our calculations are for the second output layer, at 7.56 m depth, which is representative 
of the surface motion but limits the effect of direct wind drag (we have also studied the 
layer at 97 m depth; results on this layer are briefly shown in Fig. [sj. See [26] and [25] 
for a thorough evaluation of the model performance. 

Among Lagrangian techniques used to quantify ocean transport and mixing, local 
Lyapunov methods are being widely used. The idea in them is to look at the dispersion of a 
pair of particles as they are transported by the flow. To calculate FTLEs, pairs of particles 
infinitesimally close are released and their separation after a finite time is accounted; for 
FSLEs [2j two finite distances are fixed, and the time taken by pairs of particles to separate 
from the smallest to the largest is computed. Both methods thus measure how water 
parcels are stretched by the fiow, and they also quantify pair dispersion. The methods 
can also be tailored to reveal two complementary pieces of information. On the one hand 
they provide time-scales for dispersion and stirring process [H El [9l [231 [IHl US |30] . On 
the other, they are useful to identify Lagrangian Coherent Structures (LCSs), persistent 
structures that organize the fiuid transport [T5l[l3l[8l|2ll[22l|22[Tl[371ll[nill0l[2g This 
second capability arises because the largest Lyapunov values tend to concentrate in space 
along characteristic lines which could often be identified with the manifolds (stable and 
unstable) of hyperbolic trajectories [151 [131 [m [161 137] . Since these manifolds are material 
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lines that can not be crossed by fluid elements, they strongly constrain and determine 
fluid motion, acting then as LCSs that organize ocean transport on the horizontal. Thus, 
eddies, fronts, avenues and barriers to transport, etc. can be conveniently located by 
computing spatial Lyapunov fields. We note however that more accurate characterization 
of LCSs can be done beyond Lyapunov methods [TB], that high Lyapunov values can 
correspond also to non-hyperbolic structures with high shear [12] , and that an important 
class of LCSs is associated to small, and not to large values of the Lyapunov exponents 

In the present work, however, we are more interested in obtaining the first type of 
information, i.e. in extracting characteristic dispersion time-scales, quantifying the inten- 
sity of stirring, for the different ocean regions and seasons. In particular we want to focus 
on the transport process associated to eddies and other mesoscale structures. Previous 
Lagrangian analyses of the global ocean [U |12] used FTLE to quantify such horizontal 
stirring. This quantity depends on the integration time during which the pair of parti- 
cles is followed. FTLEs generally decrease as this integration time increases, approaching 
the asymptotic value of the infinite-time Lyapunov exponent [42J. We find difficult to 
specify finite values of this integration time for which easy-to-interpret results would be 
obtained across the different ocean regions. But for the mesoscale processes on which we 
want to focus, characteristic spatial scales are related to the Rossby Deformation Radius 
(RDR), with easily defined values and latitudinal dependence (see below). Thus, we use 
in this paper FSLEs as a convenient way to identify characteristics of stirring by mesoscale 
processes. FSLE are also convenient in finite ocean basins, where relevant spatial scales 
are also clearly imposed [H, [TJ [23]. As a quantifier of horizontal stirring, measuring the 
stretching of water parcels, FSLEs give also information on the intensity of horizontal 
mixing between water masses, although a complete correspondence between stirring and 
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mixing requires the consideration of diffusivity and of the stretching directions [12] . 

More in detail, at a given point the FSLE (denoted by A in the following) is obtained 
by computing the minimal time r at which two fluid particles, one centered on the point of 
study and the other initially separated by a distance 60, reach a flnal separation distance 
6f. At position x and time t, the FSLE is given by: A(x, t, (5o, 5/) = \n{6 j / 6q) . To 
estimate the minimal time r we would need to integrate the trajectories of all the points 
around the analyzed one and select the trajectory which diverges the flrst. We can obtain 
a very good approximation of r by just considering the four trajectories defined by the 
closest neighbors of the point in the regular grid of initial conditions at which we have 
computed the FSLE; the spacing of this grid is taken equal to 60. The equations of motion 
that describe the horizontal evolution of particle trajectories are 



dt R cos 9 

d9 _ v{(f),9,t) 

~dt ~ R 



u{(f),9,t) 



(2) 



where u and v stand for the zonal and meridional components of the surface velocity 
field coming from the OFES simulations; R is the radius of the Earth (6400 km), 
is longitude and 9 latitude. Numerically we proceed by integrating Eqs. ([T]) and ^ 
using a standard, fourth-order Runge-Kutta scheme, with an integration time step dt = 6 
hours. Since information is provided just in a discrete space-time grid, spatiotemporal 
interpolation of the velocity data is required, that is performed by bilinear interpolation. 
Initial conditions for which the prescribed final separation 6f has not been reached after 
integrating all the available times in the data set are assigned a value A = 0. A possible 
way to introduce small-scale features that are not resolved by our simulated velocity fields 
is by inclusion of noise terms in the equations of motion ^ . We have recently shown [20] 
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that the main mesoscale features are maintained when this eddy-diffusivity is taken into 
account, though sub-mesoscale structures may change considerably. For global scales we 
expect the effects of noise to be even more negligible. 

The field of FSLEs thus depends on the choice of two length scales: the initial sepa- 
ration 6q (which coincides with the lattice spacing of the FSLE grid and is fixed in our 
computations to the model resolution, 6q=1/10°) and the final separation 6f. As in pre- 
vious works in middle latitudes [lOl HH [20] we will focus on transport processes arising 
from the mesoscale structures. In these studies 6f was taken about 110km, which is of 
the order of, but larger than, the mesoscale size in middle latitudes. Note that 6f should 
be a decreasing function of the latitude, since mesoscale structures decrease in size with 
Rossby Deformation Radius (RDR). We need not to exactly match RDR but to guarantee 
that our choice of 6f is similar but larger than mesoscale lengths, and also that it is a 
smooth function to avoid inducing artifacts. We have then chosen 6f as 6f = 1.3| cos^^| 
degrees; other reasonable choices lead to similar results to those presented here. 

We compute the FSLEs by backwards time integration. In this way we quantify the 
fluid deformation by past stirring. When computing LCSs this leads to structures easier 
to interpret since they can be associated with the actual shape of tracer filaments [2T|[TT]. 
However, given that forward and backward exponents in incompressible flows are related 
by temporal shifts and spatial distortions |17j, and that we are interested in temporal and 
spatial averages over relatively large scales, we do not expect significant differences when 
using forward exponents to calculate the stirring quantifiers presented below. This was 
explicitly checked in a similar framework in [TO] . 

Lagrangian measurements have been shown to correlate well with several Eulerian 
quantities at several scales 112] • In particular it is pertinent to correlate stirring with 
Eddy Kinetic Energy (EKE) since it is expected that more energetic turbulent areas 
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would also present stronger horizontal stirring, mainly due to the spawning of eddies (see 
however [Ml 133] ). Given an integration period T long enough (for instance T= one year), 
the EKE (per unit of mass) is given by: EKE = | (m'^ + v'"^), where u' and v' are the 
instant deviations in zonal and meridional velocities from the average over the period T, 
and the brackets denote average over that period. Another Eulerian measurement used 
in this work is the surface relative vorticity, given by w = |f — with positive (vs 
negative) u associated to cyclonic (vs anticyclonic) motion in the Northern Hemisphere 
(opposite signs in the Southern Hemisphere). An additional Eulerian candidate to look 
for Lagrangian correspondences is the local strain rate, but it has been shown [l3l US] to 
scale linearly with EKE^^'^ and thus it will not be explicitly considered here. 

Conditioned averages of A as a function of another variable y (let y be EKE^/^ or 



uj) introduced in Subsection 3.4 are obtained by discretizing the allowed values of y by 



binning; 100 bins were taken, each one defining a range of values (y„, ?/„+i) and represented 
by the average value tjn = ^"'''^"+-^ . So, for each discretized value of y„ the average of all 
the values of A which occur coupled with a value in (?/„, yn+i) is computed. The result is 
an estimate of the conditioned average A(?/) (which is a function of y) at the points 

3 Results 

3.1 Global horizontal stirring from FSLE 

In Fig. [T] we present a map of FSLEs at a given time. Typical values are in the order 
of 0.1 — 0.6 days^^, that correspond well to the horizontal stirring times expected at the 
mesoscale, in the range of days/weeks. Spatial structures, from filaments and mesoscale 
vortices to larger ones, are clearly identified; see a representative zoom of the South 
Atlantic Ocean (Bottom of Fig. [T|, where the typical filamental structures originated by 
the horizontal motions are evident. 
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Instantaneous maps of FSLEs have a significant signature of short-lived fast processes 
and are adequate to extract LCSs, but we are more interested in slower processes at larger 
scales. We have hence taken time averages of FSLEs over different periods, in order to 
select the low-frequency, large-scale signal. In this way we can easily characterize regions 
in the global ocean with different horizontal stirring activity; areas with larger values of 
averaged FSLEs are identified as zones with more persistent horizontal stirring [10], as 
shown in Fig. [2^. As expected, we can observe that high stirring values correspond to 
Western Boundary Currents (WBCs) and to the Antarctic Circumpolar Current, while 
the rest of the ocean and the Eastern Boundary Currents (EBCs) display significantly 
lower values. 

3.2 Geographical characterization of horizontal stirring 

A convenient quantity used to characterize stirring in a prescribed geographical area A 
was introduced by [lOj, which is simply the spatial average of the FSLEs over that area at 
a given time, denoted by < A(x, t) >a- Time series of this quantity for the whole ocean 
and the Northern and Southern hemispheres are shown in Fig. |3^. It is worth noting that 
the stirring intensity is typically larger in the Northern Hemisphere than in the Southern 
one. 

Further information can be obtained by analyzing the FSLE Probability Distribution 
Functions (PDFs). In Fig. |3]d we present the PDFs for both hemispheres and the whole 
ocean; the required histograms are built using A values computed once every week during 
one year (52 snapshots) at each point of the spatial FSLE grid in the area of interest. 
Each one of these PDFs is broad and asymmetric, with a small mode (i.e., the value of 
A at which the probability attains its maximum) and a heavy tail. Similarly to what was 
discussed by and ^5] for the FTLE case, these PDFs are well described by Weibull 
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distributions with appropriate values for the defining parameters. We note that an exphcit 
relationship between FTLE and FSLE distributions was derived by |4lj , but we have not 
checked if our flow is in the regime considered in that reference. The mode for the 
Southern Hemisphere is smaller than that of the Northern Hemisphere. Thus, Northern 
Hemisphere is globally more active in terms of horizontal dispersion than the Southern 
one. The same conclusions hold when looking at seasonally averaged instead of annually 
averaged quantities (not shown). 

Taking into account the observed differences between Northern and Southern Hemi- 
spheres, we have repeated the same analyses over the main ocean basins in a search for 
isolating the factors which could contribute to one or another observed behaviors. In Fig. 
^ we show the time evolution of < A >a as computed over the six main ocean basins 
(North Atlantic, South Atlantic, North Paciflc, South Paciflc, Indian Ocean and Southern 
Ocean), compared to the one obtained over the global ocean. The Southern Ocean hap- 
pens to be the most active (in terms of horizontal stirring) because of the presence of the 
Antarctic Circumpolar Current, followed by the Atlantic and Indian Oceans, and flnally 
the Paciflc. We have also computed (Fig. [sji) PDFs of FSLE for the different oceans. As 
before, we obtain broad, asymmetric PDFs with small modes and heavy tails. The small- 
est mode Am corresponds to the Southern Paciflc, meaning than there is less horizontal 
stirring activity in this basin, in support of what is also visually evident in Fig. [3]:. On 
the opposite regime we observe that the largest FSLE values correspond to the Southern 
Ocean. For the rest of oceans the PDFs are rather coincident with the whole ocean PDF. 

We have gone further to a smaller scale, by repeating the same analyses for the main 
currents in the global ocean: Gulf Stream, Benguela, Kuroshio, Mozambique, East Aus- 
tralian, California, Peru and Canary currents. As evidenced by Fig. [3^ there is a clear 
separation in two groups of currents in terms of their horizontal stirring properties: the 
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most active currents (including Gulf Stream, Kurosliio, Mozambique and East Australian 
currents, all of them WBCs) and the least active ones (including Benguela, California, 
Peru and Canary Currents, which correspond to EBCs). The distinction remains in the 
PDF analysis: we can clearly distinguish two groups of PDFs: a) narrow PDFs highly 
peaked around a very small value of A (EBCs); b) PDFs peaking at a slightly greater 
value of A, but significantly broader (WBCs). Since the PDFs of the WBCs are broader, 
large values of FSLEs are found more frequently, i.e., more intense stirring occurs. This 
appears to be a reflection of the well-known mechanism of Western Intensification by |39j . 
Also, the asymmetry and tails of the PDFs show that the FSLE field is inhomogeneous 
and that there are regions with very different dispersion properties. Following [3], asym- 
metry and heavy tails make the PDFs quite different from the Gaussians expected under 
more homogeneous mixing. These characteristics are then indications that chaotic motion 
plays a dominant role versus turbulent, smaller scales, dynamics. That is, the large scale 
velocity features control the dynamics, something that is also reflected in the filamentary 
patterns of the LCS shown in Fig. [T} 

3.3 Seasonal characterization of horizontal stirring 

Horizontal stirring in the global ocean has a strong seasonal variability, as shown in 
Fig. [3^. Maximum values of < A >a in the Northern Hemisphere are reached early in 
that hemisphere Summer, and minimum ones early in that hemisphere Winter. The same 
happens for the Southern hemisphere related to its Summer and Winter periods. 

Seasonally averaged FSLEs in the whole ocean over the four seasons are shown in 
Fig. |4j The spatial pattern is rather similar in all of them, and also similar to the 
annually-averaged spatial distribution shown in Fig [2^. Higher FSLE levels are found 
at the Gulf Stream and Kuroshio in the Northern Hemisphere in Spring and Summer of 
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that hemisphere. Analogously for the Eastern Australia and Mozambique Currents in the 
Southern Hemisphere relative to their own Spring and Summer time. 

Following |13], to analyze which areas are more sensitive to seasonal changes, we 
computed the standard deviation of the annual time series of FSLE (see Fig. [s]). Larger 
values appear to correspond to the more energetic regions thus showing a higher seasonal 
variability. More information about seasonal variability of different oceanic regions can 
be obtained again from Fig. |3| Time evolution of stirring in the North Atlantic and North 
Pacific, shown in Fig. |3]3, attains high values in Spring and Summer, and minimum ones 
in Winter. Concerning the main currents, we found than values of stirring in Kuroshio, 
Gulf Stream, East Australia, and Mozambique currents increase in Spring and Summer 
and decrease in Winter (see Fig. [3^). This seasonal variability is also present in EBCs 
but the amplitude of the changes is smaller than in WBCs. 

The generic increase in mesoscale stirring in Summer time detected here with Lya- 
punov methods has also been identified in previous works and several locations [HI [321 
l27l [3n m] (in most of the cases from the EKE behavior extracted from altimetric data). 
Although no consensus on a single mechanisms seems to exist (see discussion in [H]) 
enhanced baroclinic instability has been proposed in particular areas [32], [31], as well as 
reduced dissipation during Summer [S]. 

We have also calculated longitudinal (zonal) averages of the time averages of FSLE 
in Figs. [2^ and |4| This is shown in Fig. |6] (top figure for the Northern hemisphere and 
bottom figure for the Southern one). First of all, we see that horizontal stirring has a 
general tendency to increase with latitude in both hemispheres. One may wonder if this is 
a simple consequence of the decreasing value of 5/ we take when increasing latitude. We 
have checked that the same increasing tendency remains when the calculation is redone 
with a constant 5f over the whole globe (not shown), so that this trivial effect is properly 
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compensated by the factor ln(5//(5o) in the FSLE definition, and what we see in Fig. [6] is 
really a stronger stirring at higher latitudes. Note that this type of dependence is more 
similar to the equivalent sea surface slope variability, K^i, calculated from altimetry in 
[38] than to the raw zonal dependency of the EKE obtained in the same paper. Since Ksi 
is intended to represent Sea Surface Height variability with the large scale components 
filtered out, we see again that our FSLE calculation is capturing properly the mesoscale 
components of ocean stirring observed by other means. 

It is also clearly seen that latitudinal positions of local maxima of stirring correspond 
to the main currents (e.g. Gulf Stream and Kuroshio around 35°N; Mozambique, Brazil 
and East Australia around 25°S). The picture in Fig. [6] confirms that horizontal stirring 
is somehow higher in local Summer in mid-latitudes, were the main currents are, for 
both hemispheres. At low and high latitudes however the horizontal stirring is higher 
in local winter-time for both hemispheres, which is particularly visible in the Northern 
Hemisphere at high latitudes. A similar behavior was noted by [H] in the subpolar North 
Pacific and part of the subpolar North Atlantic for EKE derived from altimetry. Possible 
causes pointed there are barotropic instabilities or direct wind forcing. 

3.4 Lagrangian-Eulerian relations 

Lagrangian measures such as FSLEs provide information on the cumulative effect of fiow 
at a given point, as it integrates the time-evolution of water parcels arriving to that 
point. They are not directly related to instantaneous measurements as those provided by 
Eulerian quantities such as EKE or vorticity, unless some kind of dynamic equilibrium or 
ergodicity-type property is established so that the time-integrated effect can be related 
to the instantaneous spatial pattern (for instance, if the spatial arrangement of eddies 
at a given time gives an idea about the typical time evolution of a water parcel) or 
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their averages. EKE gives information on the turbulent component of the flow, which is 
associated to high eddy activity, while relative vorticity uj takes into account the shear 
and the rotation of the whole flow. Eventual establishment of such dynamic equilibrium 
would allow to substitute in some instances time averages along trajectories by spatial 
averages, so providing a useful tool for rapid diagnostics of sea state. Thus, we will relate 
the Lagrangian stirring (as measured by the FSLEs) with an instantaneous, Eulerian, 
state variable. Of course, the Lagrangian-Eulerian relations will be useful only if the 
same, or only a few functional relationships hold in different ocean regions. If the relation 
should be recalculated for every study zone, the predictive power is completely lost. 

We have thus explored the functional dependence of FLSEs with EKE and relative 
vorticity. In Fig. [2] the time average of these three fields is shown. Comparing FSLEs 
(Fig. |2^) and EKE (Fig. [2|d), we see that high and low values of these two quantities are 
generally localized in the same regions. There are a few exceptions, such as the North 
Pacific Subtropical Countercurrent, which despite being energetic |32] does not seem to 
produce enough pair dispersion and stretching at the scales we are considering. It was 
already shown by and [12] that variations in horizontal stirring are closely related 
to variations in mesoscale activity as measured by EKE. Note the similarity, with also 
an analogous range of values, of the EKE plot in Fig. |2}d), obtained from a numerical 
model, to that of pE] (first figure), which is obtained from altimetry data. In [43J a 
proportionality between the stretching rate (as measured by FTLE) and EKE^/^ was 
inferred for the Tasman Sea (a relation was found but no fit was attempted in the global 
data set described in [12] )• In order to verify if a similar functional dependence between 
FSLE and EKE could hold for our global scale dataset, we have computed different 
conditioned averages (see Section [2]), shown in Fig. [7j in the left panel we present the 
conditioned average X{EKE), while in the right panel A(ci;) is shown; both functions were 
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derived from the time averaged variables shown in Figure [2] 

The smooth curve depicted in Fig. [7| left, is an indication of a well-defined functional 
relationship between A and EKE, similar to the ones found by |13] and [12] from altimeter 
data. Notice however that the plot just gives conditioned averages, but the conditioned 
standard deviation -which is a measure of randomness and fluctuations- is not negligible. 
An idea of the scatter is given for selected areas in Fig. [8| Considerably less compact 
relationships were obtained in the Mediterranean sea [1^. Fig. |8] shows that very different 
dynamical regimes identified by different values of A may correspond to the same level of 
EKE. As a Lagrangian diagnostic, we believe that FSLE is more suitable to link turbulence 
properties to tracer dynamics than Eulerian quantifiers such as EKE. FSLEs provide 
complementary information since very energetic areas, with large typical velocities, do 
not necessarily correspond to high stretching regions. A paradigmatic example is a jet, 
or a shear flow, where small dispersions may be found because of the absence of chaotic 
trajectories. A functional relation between A and u is also obtained (Fig. [7| right), 
although it is much noisier and probably worse-behaved. When particularizing for the 
different regions, we see that for EKE the WBCs are all roughly associated with one 
particular functional relation for the conditioned average A while EBCs gather around a 
different one. None of the two prototype Lagrangian-Eulerian relations fits well to the 
relation A oc EKE^^^ proposed for FTLE by from altimeter data in the Tasman sea. 
Data are too scarce to make a reliable fitting for the conditioned average, in particular 
for the EBC In Fig. [8] we see that relations of the form A oc EKE"' could be reasonably 
fitted to scatter plots of the data, with a larger than the 0.25 obtained in [13], specially for 
WBC were a is in the range (0.34,0.40). This quantitative difference of our results with 
[33] may rest upon the fact that they considered just the Tasman Sea and we consider 
the different oceans. Other sources for the difference could be that we are using FSLE 
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of velocity data from a numerical model, instead of FTLE from altimetry, or that they 
use a grid of relatively low resolution 0.5° x 0.5°, while ours is 0.1° x 0.1°. Maybe their 
coarser resolution is not enough to resolve filaments which are the most relevant structures 
in our FSLE calculations. Despite this the qualitative shape of the Lagrangian-Eulerian 
relations is similar to the previous works ^\ W2\ . 

In order to analyze the ocean regions beyond boundary currents, we have also com- 
puted the conditioned averages for the Equatorial Current and for a 40° longitude by 
20° latitude sub-region centered at 245° longitude and —30° latitude in the middle of the 
sub-tropical gyre in the Pacific Ocean (and hence an area of scarce horizontal stirring 
activity). We see (Fig. [7| left) that the EBC Lagrangian-Eulerian relation is valid for 
these two areas. We have also verified that the relations derived from annually-averaged 
quantities remain the same for seasonal averages (not shown). The important point here 
is the occurrence of just two different shapes for the EKE-FSLE relations across very dif- 
ferent ocean regions, which may make useful this type of parametrization of a Lagrangian 
quantity in terms of an Eulerian one. For the relations of FSLE in terms of relative vor- 
ticity, a distinction between WBC and EBC still exists but the results are less clear and 
class separation is not as sharp as in the case of EKE (see Fig. [7| right). For instance. 
Gulf Stream and Kuroshio, despite being both WBC, do not seem to share the same 
Lagrangian-Eulerian relation, which limits its usefulness. 

4 Conclusions 

In this paper we have studied the space and time variability of horizontal stirring in the 
global ocean by means of FSLE analysis of the outputs of a numerical model. Similarly 
to what has been done in previous works, FSLEs can be taken as indicators of horizontal 
stirring. Being Lagrangian variables, they integrate the evolution of water parcels and 
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thus they are not completely local quantities. We have taken averages to analyze two main 
time scales (annual and seasonal) and three space scales (planetary scale, ocean scale and 
horizontal boundary scale). Our velocity data were obtained by using atmospheric forcing 
from NCEP. Structures and dynamics at small scales will be probably more realistic if 
forcing with higher resolution observed winds, as in |36j. But since we have not studied 
the first model layer which is directly driven by wind, and we have focused on averages at 
relatively large time and spatial scales, we do not expect much differences if using more 
detailed forcing. 

Horizontal stirring intensity tends to increase with latitude, probably as a result of 
having higher planetary vorticity and stronger wind action at high latitudes, or rather, 
as argued in [H] because of barotropic instabilities. Certainly, new studies are required 
to evaluate these hypothesis. At a planetary scale we observe a significantly different 
behavior in the Northern hemisphere with respect to the Southern Hemisphere, the first 
being on average more active in terms of horizontal stirring than the second one. This 
difference can probably be explained by the greater relative areas of subtropical gyres in 
the Southern Hemisphere with small stirring activity inside them, which compensates in 
the averages the great activity of the Antartic Circumpolar Current. At an ocean scale, 
we observe that the level of stirring activity tends to decay as the size of subtropical 
gyres increases, what is an indication that the most intense horizontal stirring takes place 
at the geographical boundaries of ocean basins. For that reason, we have finally ana- 
lyzed the behavior of stirring at boundary scale, which is mainly related to WBCs and 
EBCs. EBCs behave in a similar way to ocean interior in terms of all the quantities we 
have computed, including the Lagrangian-Eulerian relations. Thus, the main hot spots of 
horizontal stirring in the ocean are WBC. The observed small mode in the global FSLE 
PDFs also indicates that horizontal stirring is not very intense for the vast majority of 
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the ocean, but the heavy tails indicate the existence of large excursions at some specific, 
stretched locations (e.g., inside the WBCs and other smaller scale currents active enough 
to generate stirring). This type of uneven distribution is characteristic of multifractal sys- 
tems arising from large scale chaotic advection, something that was discussed for oceanic 
FSLEs in [20] • 

Regarding seasonal variability, generally we observe stronger stirring during each hemi- 
sphere's Summer time. Medium and high latitudes behave however in the opposite way: 
stirring is more active during the hemisphere Summer for medium latitudes and during 
the hemisphere Winter for high latitudes. Medium latitudes are strongly affected by the 
behavior of WBC, which experience intensification of horizontal stirring during Summer 
[TS| EZl ED Hi]- As commented before, high latitude Winter intense stirring could be 
the result of a stronger action of wind during that period or of barotropic instabilities 
[S], and dedicated studies are required to evaluate these hypothesis. 

Finally, we have studied the connection between time-extended Lagrangian FSLEs 
and instant Eulerian quantities such as EKE and relative vorticity. For the case of EKE, 
the different ocean regions give rise to just two different Lagrangian-Eulerian relations, 
associated to an intense or a weak stirring regimes. The existence of these two regimes 
implies that pair dispersion and stretching strength are larger in a class of ocean areas 
(represented by WBC) than in another (e.g. EBC) at mesoscales, even when having the 
same EKE. 
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Figure 1: Top: Snapshot of spatial distributions of FSLEs backward in time corresponding 
to November 11, 1990 of the OFES output. Resolution is 6q = 1/10°. Bottom: Zoom in 
the area of the box inside top figure (South Atlantic Ocean). Coherent structures and 
vortices can be clearly seen. The colorbar has units of day~^. 
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Figure 2: a) Time average of the FSLEs in the Global Ocean. Geographical regions of 
different stirring activity appear. The colorbar has units of day~^. b) Spatial distribution 
of annual EKE^^'^ (cm/s). c) Time average of Relative Vorticity (w) in the Global Ocean. 
The color bar has units of day~^ . In all the plots the averages are over the 52 weekly 
maps computed from November 1st, 1990 to October 31th, 1991. 
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Figure 3: Left column: Temporal evolution (from November 1st, 1990 to October 31th, 
1991) of the horizontal stirring (Spatial average of FSLEs). Right column: PDFs of the 
FSLEs (histograms are built from the A values contained at all locations of the 52 weekly 
maps computed for the second simulation output year). Top: for both hemispheres and 
for the whole ocean. Middle: for different oceanic regions. Bottom: for some main 
currents during one simulation year. In addition to the results from the second surface 
layer analyzed through the paper, panel a) shows also stirring intensity in a layer close to 
100m depth. 
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Figure 4: Time average of the FSLEs in the Global Ocean for the each season. Spring: 
from March 22 to June 22. Summer: from June 22 to September 22. Autumn: from 
September 22 to December 22. Winter: from December 22 to March 22. The colorbar 
has units of day~^. 
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Figure 5: Standard deviation of weekly FSLE maps of one year. The colorbar has units 
of day~^ 
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Figure 6: Cross-ocean zonal average of the annual, relative Summer and relative Winter 
time average of FSLE maps from Fig|2^ as a function of latitude (expressed as absolute de- 
grees from Equator to make both hemispheres comparable). Top: Northern Hemisphere; 
bottom: Southern Hemisphere. 
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Figure 7: Left: Lagrangian-Eulerian relations. Left: the conditional average A e 
function of its corresponding annually averaged (second year) EKE for different regions 
and currents. We clearly observe two groups of relations FSLE-EKE. Right: same plot for 
the conditional average X^j as a function of its corresponding annually averaged (second 
year) u. Although we observe also the same two two groups of FSLE-o; relations, these 
functions are much noisier and region- dependent. 
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Figure 8: Scatter plots showing temporally averaged FSLE values at different spatial 
points in regions of Fig. [2^, and EKE values (as displayed in Fig. ^p) at the same points. 
The regions displayed here are eight of the main currents. Fittings of the type y = cX^ 
are also displayed, where y is the temporal mean of FSLE and X is EKE^^^. Note that 
this implies < FSLE >= c EKE" with a = b/2. 
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